Dynamical Casimir Effect in cavity with AMevel detector or iV — 1 two-level atoms 
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We study the photon generation from vacuum via the Dynamical Casimir Effect in a cavity 
containing a iV-level detector in equally-spaced resonant ladder configuration or N — 1 identical 
resonant two-level atoms. If the modulation frequency equals exactly twice the unperturbed cavity 
frequency, the photon growth goes on steadily for odd N, while for even N at most N — 2 photons 
can be generated. This finding is corroborated by numerical calculations. In the limit TV — > oo we 
obtain the harmonic oscillator model of the detector, which admits analytical results. 
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PACS numbers: 42.50.Pq, 32.80.-t, 42.50.Ct, 42.50.Hz 



A possibility of creating quanta of the electromagnetic 
field from the initial vacuum state in cavities with moving 
boundaries, called nowadays as the Dynamical Casimir 
Effect (DCE), was a subject of numerous theoretical 
studies for a long time: see, e.g., the most recent reviews 
PHI]. It was shown Q that one might expect a signifi- 
cant rate of photons generation inside ideal cavities with 
resonantly oscillating boundaries. The simplest model 
describing this effect takes into account a single reso- 
nant cavity mode whose frequency is rapidly modulated 
according to the harmonical law uj t = u>o[l + esin(^t)] 
with a small modulation depth, |e| -C 1. We shall use 
dimensionless variables, setting K = ojq = 1. Then the 
Hamiltonian for the resonance mode has the form Q 



H c =w t n- ixt(a 2 
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), xt = (4ut)- l dut/dt, (1) 



where a and are the cavity annihilation and creation 
operators, and n = a^a is the photon number opera- 
tor. It is well known that the number of photons created 
from the initial vacuum state is maximal if the modu- 
lation frequency is exactly twice the unperturbed mode 
frequency, i.e., rj = 2. The mean number of photons (n) 
and the Mandel factor Q = [((An) 2 ) — (n)]/(n) increase 
with time in this ideal case as (hereafter we use the sub- 
script for the quantities related to the empty cavity) 

(no(i)) = sinh 2 (rf/2), Qo(t) = 1 + 2(n {t)). (2) 

The field mode goes to the squeezed vacuum state with 
the following variances of the field quadrature operators 
x = (a + at)/v2 and p = (a — aJ)/(y/2i) (the average 
values of these operators are zero) 
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But simple formulas © and ([3]) hold for the ideal 
empty cavity only. To registrate the emerging photons 
one has to couple the field mode to some detector. And 
here the problem of the back action of the detector on the 
field arises, because in many realistic cases the coupling 
between the field and detector can be much stronger than 
that between the field and vibrating cavity walls. This 
was noticed in Q, where it was shown that for the sim- 
plest model of detector as a two-level "atom" , no photons 



can be created at all for the modulation frequency rj = 2, 
if the field atom coupling constant g is much bigger than 
the frequency modulation amplitude s. A more detailed 
investigation of this problem was given recently in [3, H| , 
where different resonant regimes were found and analyzed 
for different ratios gje (both small and big). On the other 
hand, in papers [3, [10| we have discovered that when the 
inter-level transitions are resonant with the cavity bare 
frequency and the modulation frequency equals twice the 
cavity unperturbed frequency, many photons can be gen- 
erated for the 3-level atom in ladder configuration or a 
chain of two 2-level atoms, even if g ^> e. Therefore, 
a question naturally arises: how many photons can be 
generated if one uses the AMevel atom in the equidis- 
tant ladder configuration and the inter-level transitions 
are resonant with the cavity mode modulated at exactly 
twice the unperturbed cavity eigenfrequency? We answer 
this question in this brief report. 

If the selected field mode interacts with the TV-level de- 
tector in resonant ladder configuration, the Hamiltonian 
describing the whole system "field mode + atom" can be 
taken in the Rotating Wave Approximation (RWA) form 



N 



JV-1 

E 



H = H c + 22 E i a H + Z2 ffi( a<T »+l,« + at(7 M+l)> ( 4 ) 



where Ei is the energy of the i-th atomic eigenstate |i) (in 
bold), Gi.j = |i)(j| is the generalized Pauli operator, and 
gi (assumed to be real) is the coupling parameter between 
the atomic states {|i), |i + 1)} through the cavity field. 
The condition of validity of RWA is \gj\ <C 1. We consider 
here the resonant case, so we assume -Ej+i — Ei = 1 for 
i = l,..., N — 1. Our aim is to find out how the photon 
generation from vacuum is affected by the detector in the 
weak modulation regime, |e| <C gj (j = 1, . . . , N — 1) for 
the modulation frequency rj = 2. 

To find the wave function of the whole system \^(t)), 
we make the transformation ^(t)) = V(t)\ip(t)) with 
V(t) = exp[-it(n+Y^ =1 E l a ll )]. Then after RWA we get 
the following Hamiltonian governing the time evolution 
of \tp(t)) in the interaction picture (where j3 = e/4): 
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Hi = -i/3 (a 2 - a t2 ) + ^ gi(aa i+ i ti + aV iii+ i). (5) 
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In the special case when N —> oo and gj = \fjg the 
detector becomes a simple harmonic oscillator (HO) if we 
associate the atomic level |1) with the oscillator ground 
(zero energy) state and make the replacements 

oo oo 
3=1 3 = 1 

where b and are the annihilation and creation opera- 
tors associated with the detector ([&, b^] = 1). Then the 
Heisenberg equations of motion for operators a(t) and 
b(t) can be solved exactly: 

a(t) = (Cpc 7 + pSps^) a>o — igCpSjbo 

+ (Spay + ftCps^al+igSps^bl, (6) 



(n(t)), Q(t) and (Ap) 2 (t), which are clearly seen in Figs. 
CD and [U below. For fit > 1 Eq. © can be simpli- 
fied as Q(t) ~ 2(n(t)), which is a typical relation for 
highly squeezed vacuum states. Nonetheless, the state of 
the field mode is not exactly the vacuum squeezed one, 
since the uncertainty product A = (Ap) 2 (Ax) 2 is bigger 
than the minimal possible value 1/4: A = 1/4 + g 2 /3 2 s 7 
(the covariance between x and p quadratures is zero in 
the case discussed). As soon as the initial state of the 
system was Gaussian, it remains Gaussian for all times. 
In such a case, the purity of the field mode equals [l3[ 
fx = Tr(/5 2 ) = (4A) -1 / 2 , where p is the statistical opera- 
tor. If \g\ ^> |/3 1, the purity is only slightly below unity, 
oscillating with the frequency 2j. But for < \g\ < |/?| 
the purity goes monotonously to zero. 



b(t) 



{CpCj - PSpSy) b - igCps 7 a 

(SpCy - PCpSj) bl - igSps^al, (7) 



where 7 = \J g 2 — (3 2 , Cp = cosh(/3i), Sp = sinh(/3i), 
c 7 = cos(7t), and s 7 = sin(7t)/7. If \(3\ > \g\, then 
trigonometrical functions should be replaced by their hy- 
perbolic counterparts with 7 replaced by 7 = \//3 2 — g 2 . 
The mean number of quanta in the field mode (n(t)} and 
the mean excitation number of the detector (rib(t)) for 
the initial vacuum state are equal to 
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The Mandel factor and quadrature variances of the field 
mode are 
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lfp = 0, thenEqs. ©-(21 go to ©-©. Butif| 3 |> |/3|, 
then the rate of photon generation becomes roughly twice 
smaller than in the empty cavity (16J. Moreover, the ratio 
(nb(t)) / (n(t)) is close to unity if |^| » \(3\. For e > and 
/3t 3> 1 we see the exponential growth of the mean photon 
number with increment e/2, which is modulated by some 
oscillations with the frequency 27: 



(n(t)) 



1 



,2/3* 



1 



■7 



sin(27t) 



2^ 

7 2 



sin 2 (7t) 



(11) 



For t = t n + 5t with 2jt n = (2n + 1)tt we have the 
following Taylor expansion of formula (jlip : 

(n(t„ + 5t)) « (n(t n )) [1 + (4/3)/3 3 2 (^) 3 ] . 

Consequently, the function (n(t)) is practically constant 
in some neighborhood of i„. This results in the ap- 
pearance of almost horizontal "shelves" in the plots of 




FIG. 1: (Color online) Various quantities as functions of the 
dimensionless time et for different numbers of the detector 
levels for gj = \fjg with g = 10 -2 and e = 10 -3 . Blue thick 
lines represent the results for harmonic oscillator detector, 
while the red dashed line - for the empty cavity. 



For equidistant detectors with a finite number of levels 
analytical solutions cannot be obtained. In this case we 
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solve the Schrodinger equation for the Hamiltonian ([5]) by 
expanding the wavefunction in the Fock basis as \ip(t)) = 

YljLi YlkLaPj,k(t)\ji k), where the first index stands for 
the atomic eigenstate and the second for the Fock state 
of the cavity field. We get the set of differential equations 



Pj,m = PW m(m - l)p Jim _ 2 - y/(m + l)(m + 2)p^ m+2 ) 
-i{gjV^Pj+i,m-i + 9j-iVm~+Tpj-i,m+i) (12) 

where go — gj>N = and the dot stands for the time 
derivative. Hereafter we suppose that |/3| <C \gj\- To 
understand qualitatively the resulting dynamics in this 
case we follow the method employed in 0, Gil • First 
we solve the equations (TT2"j) for (3 = 0, obtaining solutions 
in the form of exponentials with time- varying arguments 
multiplied by constant coefficients. Then we substitute 
the solutions obtained back into (fl2|) with (3 ^ 0, as- 
suming that now the coefficients at the exponentials be- 
come time-dependent. Thus we obtain a set of differential 
equations for these coefficients. 

In this paper we assume that the only nonzero initial 
probability amplitude is pi^o(0) = 1, so the first equation 
to be solved is Pi,o(t) = — V^PPi,2(t)- Our first task is 
to find the expression for pi i2 (t) when j3 = 0; it must be 
obtained from the matrix equation 



X = -MX , 

where X = (pi, 2 ;P2,i;P3,o) and 

y/2gx 
M = ( V2 5l g 2 
52 



(13) 



(14) 



For N = 2 the eigenvalues of matrix M are ±v2£fi, 



so Pi, 2 = J2k=i ^k' ex P{-iV K k' t ), where A [ ^' are con- 
stant coefficients and ip k are the eigenvalues of M. 
Since these eigenvalues are big, \<p) J | \/3\, no sig- 

(2) 

nificant coupling between the coefficients p\^ and A, 
(with k — 1,2) arises when one puts the ansatz pi j2 — 

J~f k=1 A (2 \t)exp(-iip (2) t) into Eq. (I2J with /3 ^ 0. So 
within the RWA reasoning one has pi^o — 0. Thus p\^ 
is decoupled from all the other probability amplitudes 
and the system remains in the initial state |1, 0). On the 
other hand, for N > 2 the solution of (fT3"|) for p lj2 is 

Pi.2 = Ylk=i A-k cx p(~ *Vk w ith the third eigenvalue 

(3) 

of matrix M equal to zero, <y9 3 = 0. Moreover, one can 

(3) 

check that p 2 ,i does not contain the term A3 , but p3 t o 
does (exact expressions are given in Q). Hence, in this 

(3) 

case pi : o is resonantly coupled to A 3 , so pi j2 and 753,0 
become populated, while p 2 .\ is coupled off-resonantly, 
so it is much smaller than the former two coefficients 
(|A^ 3) [ 5 [A£ 3) ] < |A| 3) |). Thus in this case at least two 
photons can be created for sure. 

,(3) 



(2), 
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I3J with X = (pi,4;P2,3;P3,2;P4,i;P5,o) and 

/ 51V4 \ 
giV^ g 2 y/3 

M= 52^3 53V2 
53^2 54 

V 54 j 



For N — Awe have pi.4 = Ylt=i ex P(~ ^fc 4 ^)' where 
all the eigenvalues are different from zero and much larger 

(3) 

than being functions of (i — 1, 2, 3). Therefore A 3 

does not couple resonantly to any A^ . Thus at most two 
photons are generated for the four-level atom. For N > 
4 the solution for p\^ is p\^ — J2k=i ex P( — Wk^) 



( 5^ in 

with ip 5 — 0, and A^ also appears in p% 2 and p$ 0, 

(3) ^ 



but not in p 23 and P41. Thus, the coefficient A 

(5) 

resonantly coupled to A 5 and at least four photons can 
be generated. 

Continuing this analysis for higher order matrices, one 
can easily verify that matrices M of odd orders always 
have one null eigenvalue, and the coefficient A that multi- 
plies the exponential with the null eigenvalue appears in 
probability amplitudes with odd atomic level and even 
photon number. On the other hand, all eigenvalues of 
matrices M of even orders are different from zero and 
much larger than \/3\. Hence we arrive at the general 
rule: The number of created photons is unlimited for odd 
numbers of levels, while at most N — 2 photons can be 
created if the number of levels N is even. Moreover, 
the probability of detecting an odd number of photons 
is much smaller than the probability of detecting an even 
number of photons. 

These results can be extended straightforwardly to an- 
other practical scenario: (N — 1) resonant 2-level atoms 
interacting with the field mode, when initially all the 
atoms are in their ground states and the field is in the 
vacuum state, for the modulation frequency 77 = 2. It was 
shown in [loj that for a single atom no photons are gen- 
erated for this modulation frequency, and for two atoms 
the photons are generated steadily provided the coupling 
strengths are equal. So here we assume that all the 
couplings are equal to 5. In the interaction picture the 
Hamiltonian describing this setup is the Tavis-Cummings 
Hamiltonian [10t Il4j with additional parametric amplifi- 
cation term analogous to Eq. ([5]) 



Hi = -i(3{a 2 - a) 2 ) + g(aS + + o + 5_), 



(15) 



where S± = YljLi^ °± are the collective ladder op- 
erators and aj = |lj-)(2j|, cr+ = |2 3 -)(lj| are the 
standard 2-level Pauli operators, where and \2f) 
are the ground and excited states of the j-th atom 
(j = 1,...,N — 1), respectively. Writing the wave- 
function associated to the Hamiltonian (IT5|) as \ip(t)) = 

N 



Now we must see the coupling of the coefficients A^ ' Ej=i EfcLoft.felj' where |j) denotes the symmetric 
(j = 1, 2, 3) to the next subset, so we have to solve Eq. normalized Dicke state [l5| with (j — 1) excitations (j = 
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1, . . . ,N) [13] and \k) is the cavity Fock state, and us- 
ing the known properties S+|j) = \fj (N — j)\j + 1) and 
S-\j) = -1)(N -j + l)|j - 1), we obtain precisely 
the equation (JT2J) for the probability amplitudes, where 
the effective coupling strengths are gj = gy/ j {N — j) 
(j = 1,,..,N). Thus our statement holds unaltered. 
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FIG. 2: (Color online) Same as Fig. [TJfor g, = y/j (N — j)g. 

We confirmed these results by performing numerical 
simulations in which we solved the equations (fT2"j) us- 
ing the Runge-Kutta-Verner fifth-order and sixth-order 
method for different forms of the coupling strengths [l8[ . 
In Fig. [T]we show the behavior of the mean photon num- 
ber in), the Mandel factor Q, the squeezed quadrature 
variance (Ap) 2 and the probability P\ of finding the de- 
tector in the ground state as function of dimensionless 
time et for different number of detector's levels N. We 
verified that the probability of finding the detector in an 
even energy level is always very low (data not shown). In 
Fig. [TJwe set gj = y/Jg (j = 1, . . . ,N — 1), as one expects 
that for N ^> 1 the results should approach the behavior 
of the case of harmonic oscillator detector (depicted with 




FIG. 3: (Color online) Probability Pk of detecting K photons 
as function of dimensionless time et for N = 12, g — 10~ 2 and 
e = 1(T 3 . a) 9j = b) gj = y/j (N - jjg. 



thick blue line), and the red dashed lines denote the re- 
sults for the empty cavity according to Eqs. In 
Fig. [2] we set gj = y/j (N — j)g corresponding to N — 1 
two-level atoms. In Fig. [3] we plot the probabilities Pk 
of detecting K photons as function of time for gj = \fjg 
(Fig. W) and g, = y/j (N - Jjg (Fig. 0b) for N = 12, 
to show that at most 10 photons are created, whereas 
the probabilities of detecting odd photon numbers are 
very small (the probabilities with K > 12 are below the 
horizontal axes). Hence the numerical results are in full 
agreement with our general rule. 

In conclusion, we studied the photon pair creation from 
vacuum via the Dynamical Casimir Effect in a cavity 
containing a resonant JV- level detector in equally-spaced 
ladder configuration or N — 1 identical resonant two- level 
atoms. We considered the regime of weak modulation, 
where the cavity modulation depth is much smaller than 
the atom-field coupling. We showed that for an odd num- 
ber of levels pairs of photons are generated steadily as 
time goes on, while for an even number of levels at most 
N — 2 photons are created. For N — > oo the difference 
between even and odd values of N eventually disappears, 
and for the special choice of the coupling constants we ar- 
rived at the harmonic oscillator model of the detector. In 
this special case the rate of photon generation is roughly 
twice smaller than for an empty cavity. These findings 
can be useful for constructing efficient detectors for mon- 
itoring the photon generation via Dynamical Casimir Ef- 
fect, because the detector can be excited without inter- 
rupting the photon growth. 
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